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Abstract. In this paper we present the first results on substructuring methods for nonlocal 
operators, specifically, an instance of the nonlocal p-Laplace operator. We present a nonlocal vari- 
ational formulation of this operator, proving a nonlocal Poincarc inequality and upper bound to 
establish a spectral equivalence. We then introduce a nonlocal two-domain variational formulation 
utilizing nonlocal transmission conditions, and prove equivalence with the single-domain formula- 
tion. A nonlocal Schur complement is introduced. We establish condition number bounds for the 
nonlocal stiffness and Schur complement matrices. Supporting numerical experiments demonstrating 
the conditioning of the nonlocal single- and two-domain problems are presented. 
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1. Introduction. Domain decomposition methods where the subdomains do 
not overlap are called substructuring methods, reflecting their origins and long use 
within the structural analysis community [29]. These methods solve for unknowns 
only along the interface between subdomains, thus decoupling these domains from 
each other and allowing each subdomain to then be solved independently. One may 
solve for the primal field variable on the interface, generating a Dirichlet boundary 
value problem on each subdomain (these are Schur complement methods, see [39] 
and references cited therein), or solve for the dual field variable on the interface, 
generating a Neumann boundary value problem on each subdomain (these are dual 
Schur complement methods, see [21, 22, 20, 31]). Hybrid dual-primal methods have 
also been developed [26]. 

As domain decomposition methods are frequently employed on massively parallel 
computers, only scalable methods are of interest, meaning that the condition number 
of the interface problem does not grow (or, only grows weakly) with the number of 
subdomains. Scalable or weakly scalable methods are generated by application of an 
appropriate preconditioner to the interface problem. This preconditioncr requires the 
solution of a coarse problem to propagate error globally; see any of the references 
[6, 7, 8, 9, 19, 28, 27, 15]. For a general overview of domain decomposition, the reader 
is directed to the excellent texts [39, 30, 40]. 

All of the methods referenced above have in common that they are domain de- 
composition approaches for local problems. In this paper, we propose and study a 
domain decomposition method for the nonlocal Dirichlet boundary value problem 

£(u)=b(x), x£ft, (1.1) 
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where 

C(u):= f C(x,x')(u(x')-u(x))dx', (1.2) 

£1 is a bounded domain, Bfl is given in (2.1), b is given, and u(x) £ M. d is prescribed 
for x £ R d \f2. We prescribe the value of u(x) outside fl and not just on the boundary 
of 0, owing to the nonlocal nature of the problem. 

Nonlocal models arc useful where classical (local) models cease to be predictive. 
Examples include porous media flow [13, 14, 38], turbulence [5], fracture of solids, 
stress fields at dislocation cores and cracks tips, singularities present at the point of 
application of concentrated loads (forces, couples, heat, etc.), failure in the predic- 
tion of short wavelength behavior of elastic waves, microscale heat transfer, and fluid 
flow in microscale channels [18]. These are also cases where microscale fields are nons- 
mooth. Consequently, nonlocal models are also useful for multiscale modeling. Recent 
examples of nonlocal multiscale modeling include the upscaling of molecular dynam- 
ics to nonlocal continuum mechanics [32], and development of a rigorous multiscale 
method for the analysis of fiber-reinforced composites capable of resolving dynamics 
at structural length scales as well as the length scales of the reinforcing fibers [1]. 
Progress towards a nonlocal calculus is reported in [25]. Development and analysis 
of a nonlocal diffusion equation is reported in [2, 3, 4]. Theoretical developments for 
general class of integro-diffcrcntial equation related to the fractional Laplacian are 
presented in [10, 11, 37]. We discuss in §2 some specific contexts where the nonlocal 
operator C appears, and the assumptions placed upon C by those interpretations. 

To the best of authors' knowledge, this paper represents the first work on domain 
decomposition methods for nonlocal models. Our aim is to generalize iterative sub- 
structuring methods to a nonlocal setting and characterize the impact of nonlocality 
upon the scalability of these methods. To begin our analysis, we first develop a weak 
form for (1.1) in §3. The main theoretical construction for conditioning is in §4. We 
establish spectral equivalences to bound the condition numbers of the stiffness and 
Schur complement matrices. For that, we prove a nonlocal Poincare inequality for the 
lower bound and a dimension dependent estimate for the upper bound. We also reveal 
a striking difference between the local and nonlocal problems. Namely, the condition 
number of the discrete nonlocal operator is mesh size independent. Moreover, unlike 
the local case, we show that the condition number varies with the spatial dimension. 
In §5 wc formulate two-domain strong and weak forms of (1.1), giving particular at- 
tention to the nonlocal transmission conditions. We then prove equivalence of the 
one-domain and two-domain strong forms, and also equivalence of the one-domain 
and two-domain weak forms. Another remarkable result is that this equivalence of 
strong and weak forms can be achieved while allowing u to be less regular than is 
required to demonstrate equivalence of two- and one-domain problems for the (local) 
Laplacian. In §6, we first define an energy minimizing extension, a nonlocal analog of 
harmonic extension in the local case, to study the conditioning of the Schur comple- 
ment in the nonlocal setting. We discretize our two-domain weak form to arrive at a 
nonlocal Schur complement. We perform numerical studies to validate our theoretical 
results. Finally in §7, we draw conclusions about conditioning and suggest future 
research directions for nonlocal domain decomposition methods. 

2. Interpretations of the Operator C. The operator C appears in many dif- 
ferent application areas, from evolution equations for species population densities [12] 
to image processing [24]. We review two specific contexts in which the operator C 
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Fig. 2.1. Typical domain for (1.1). u is prescribed in BQ, and we solve for u in Q. 

of (1.1) is utilized, paying special attention to associated assumptions these interpre- 
tations place upon C in C. In all cases, we find C to have local support about x, 
meaning that we must prescribe Dirichlet boundary conditions only for 

Bil := supp(C)\ft, (2.1) 

as depicted in Figure 2.1. 

2.1. Nonlocal Diffusion Processes. The equation 

u t {x,t) =£(u(x,t)) (2.2) 

is an instance of a nonlocal p-Laplace equation for p = 2, and has been used to 
model nonlocal diffusion processes, see [4] and the references cited therein. In this 
setting, u(x, t) £ K is the density at the point x at time t of some material, and 
C(x, x') = C(x — x') is the probability distribution of material movement from x' to 
x. Then, J Rd C(x' — x)u(x', t)dx.' is the rate at which material is arriving at x from all 
other points in supp(C), and — J Rd C(x' — x)it(x', t)dx.' is the rate at which material 
departs x for all other points in supp(C) [23, 4]. 

In this interpretation of (1.1) the following restrictions are placed upon C in C. 
It is assumed that C : R d — > R is a nonnegative, radial, continuous function that is 
strictly positive in a ball of radius 6 about x and zero elsewhere. Additionally, it is 
assumed that f n C(£)d£ < oo. 

2.2. Nonlocal Solid Mechanics. The equation 

u tt (x,t) = £(u(x,t))+b(x) (2.3) 

is the linearized pcridynamic equation [34, eqn. (56)]. The corresponding time- 
independent ( "peristatic" ) equilibrium equation is (1.1). Pcridynamics is a nonlocal 
reformulation of continuum mechanics that is oriented toward deformations with dis- 
continuities, see [34, 35, 33] and the references therein. In this context, u £ R d is the 
displacement field for the body f2, and C(x, x') is a stiffness tensor. 

In this interpretation of (1.1) the following restrictions are placed upon C in C. 
It is assumed that C is bounded and strictly positive definite in the neighborhood of 
x, 7i x , defined by the following with 5 > 0: 

H x := {x : ||x'-x|| < 6}. (2.4) 

These assumptions are made to ensure material stability [34, pp. 191-194]. It is also 
assumed that C = for ||x' — x|| > 5. If the material is elastic, it follows that C(x, x') 
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is symmetric (e.g., C(x, x') T = C(x,x')). Further, it is assumed that C is symmetric 
with respect to its arguments (e.g., C(x, x') = C(x',x)). This follows from imposing 
that the integrand of (1.2) must be antisymmetric in its arguments, e.g., 

C(x,x') • [u(x') - u(x)] = -C(x',x) • [u(x) - u(x')] 

in accordance with Newton's third law. 

3. A Nonlocal Variational Formulation. Here we present a variational for- 
mulation of the nonlocal equation (1.1). For peridynamics, this was presented by 
Emmerich and Weckner in [17]. An analogous expression also appears in [34, eqn. 
(75)], as well as [25]. 

We will utilize the function space 

V := 4 (fi) = {v e Li(Sl) : v| Bn = 0} , (3.1) 

and the inner product 

(u,v)n := / u • v dx. 
J n 

The weak formulation of (1.1) is the following: Given b(x) E V, find u(x) £ V such 
that 

o(u,v) = (b,v) Vve^, (3.2) 

where 

a(u,v) :=-[[( C(x,x') • [u(x') - u(x)] dx'} • v(x) dx. (3.3) 
Jn UauBn J 

It follows from assumptions in §2 that the iterated integral in (3.3) is finite: 

- f ( f C(x,x') • [u(x') - u(x)] dx' 1 ■ v(x) dx < co, 
Jn Unuen J 

and that C(x, x') • [u(x') — u(x)] is anti-symmetric in its arguments. Combining these 
observations with Fubini's Theorem gives the identity 

- [ { [ C(x,x') • [u(x') - u(x)] dx' 1 • v(x) dx = (3.4) 
Jn UnuBn J 

\ f { [ C(x,x') • [u(x') - u(x)] dx' 1 • [v(x') - v(x)] dx. 
z Jn UnuBn ) 

Remark 3.1. The equivalence in (3.4) can be interpreted as a duality pairing. 
The equivalent expression of a(u, u) from (3.4), 



a(u, u) = - / / C(x, x') • [u(x') - u(x)] 2 dx' dx, (3.5) 
2 Jn JnuBn 



is positive. This immediately establishes coercivity. In the setting of nonlocal solid 
mechanics, (3.5) has the natural interpretation of the energy stored in deformed elastic 
solid. Note that the zero energy modes, or rigid body modes, are excluded by enforcing 
Dirichlct boundary conditions. 
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Fig. 3.1. Domain of integration for ID problem where Q = [0, 1]. The integrand is zero outside 
the region bounded by the dotted line, and is discontinuous across this line. Unlike x, the variable 
x' is integrated over BQ = [— <5, 0] U [1, 1 + 8] where u is prescribed as a nonlocal Dirichlet boundary 
condition. 



In ID with £1 := [0, 1], the weak form (3.5) becomes 

^ rl rx+5 

a(u,u) = — / / C(x, x') (u(x') — u(x)) dx'dx, (3-6) 
2 Jo Jx-S 

where the limits of integration have been adjusted to account for the support of 
C(x,x'), which is assumed to vanish if ||x — x'\\ > 8. For this problem, the two- 
dimensional domain of integration is the parallelogram shown in Figure 3.1. For 2D 
and 3D problems, the domains of integration are four and six dimensional, respec- 
tively. 

4. Nonlocal Condition Number Estimates. In the local setting, the classical 
condition number estimates rely on a Poincare inequality and an inverse inequality 
for the lower and upper bound, respectively. Similarly to the local case, we develop 
a nonlocal Poincare inequality to be used in the lower bound. We prove a Poincare 
inequality involving an explicit (5-quantification, which is a more refined result than 
what is available in the literature [4]. The (5-quantification is an essential feature in 
the nonlocal setting because the lower bound turns out to be dimension dependent, 
unlike in the local case. This dimensional dependence is induced by the neighborhood 
Tt x (see (2.4)), which is d-dimcnsional in the nonlocal setting but zero-dimensional 
(a point) in the local setting. Dimension dependence in the Poincare inequality is 
captured by S m (see §4.1) where the power m exhibits a dimensional dependence (i.e., 
m = m(d)). 

For the upper bound, we prove a direct estimate instead of an inverse inequality. 
Neither the upper bound estimate nor the Poincare requires discretization. Hence, 
our estimate is valid in infinite dimensional function spaces, a stronger result than 
that for the local setting. 

Furthermore, as in the classical condition number analysis, the dependence on 
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the kernel function C is absorbed into the following constants: 

C:= sup ||C(x,x')|| 

xgf!,x'gf!USn 

C:= inf ||C(x,x')||. 

xen.x'enueo 

Therefore, we reduce the analysis to the canonical kernel function x<5( x , x ') whose only 
role is the representation of the neighborhood in (2.4) by a characteristic function. 
Namely, 



Wxx')-| " x - x 'll^ 



(4.1) 



For the remainder of this paper, will restrict our discussion to scalar problems, e.g., 
u(x) = u(x), C(x, x') = C(x, x'), etc. 

4.1. Nonlocal Poincare inequality. We prove a nonlocal Poincare inequality 
that is valid for w(x) S L p (fl) with p > 1. Wc define a generalized bilinear form to 
support the proof in L p (Q): 

a^(u g ,u) := [ f X«(x,x') |u s (x') - u(x)|* dx'dx. 
Jn Jnusn 

For the explicit 5-quantification, we utilize inscribing and circumscribing spheres of 
il so that the Poincare constant Ap rlcr remains independent of S; see Appendix A. 
Proposition 4.1. Let fl C M. d be a bounded, simply connected domain. Assume 

that u 6 ip(^) and <? £ L p (Bfl) with m 9 (x) := < } /' „„ with p > 1. Then, 



#(x), xeSSl 

/or sufficiently small 5, there exists Xp n cr = Ap ncr .(m,p, O) > u>i</i m > 1 suc/i £/iai 

XpncrS™ I H*)\ p d*<a^\u g ,u) + 5 m ( | 5 (x')| p dx'. (4.2) 

Jen 

Proof. We construct a finite covering for f2 using strips of width ^- in a similar 
fashion utilized in [4, Prop. 2.5]: 

B :={xe£f>: ||(x,fi)|| < (4.3) 

Bi:={xen:||(x,iJo)||<^}, (4.4) 

B,:={xef2\|J :\\(^Bj-{)\\<— }, j = 2,3,...,l. (4.5) 
fc=i 

Here, for a fixed <5, the number of strips I = l(n,m). We trivially have the following 
for j = 1, 

[ [ X5(x,x , )K(x')-u(x)p' dx'dx > f f X5(x,x')|w g (x')-u(x)|f dx'dx. 
JnJnuen Jb^Jb,^ 

Note that for x' G -Bj-i and xeBj, due to £ < 1, we have 

llx-x'11 < (5 m < 5. 
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Wc refine the result in [4] by establishing an explicit 5-quantification. For that, 
we utilize the inscribing and circumscribing spheres of £1 whose radius are denoted 
by radius ri n and r ou t, respectively. Then, we compare the volume of the strips Bj 
to the volume of the annuli of width — of the inscribing and circumscribing spheres. 
For sufficiently small 5, we prove the following in Appendix A: 

c tn S m < \Bj\ <c out S m , (4.6) 

d onl 

From (4.1) we have: 



where Ci„ and c ou t depend only on f2. In fact, Ci„ ~ rf n 1 and c ou t ~ r^/ 



X s{x,x') = l, x' /.',, ;. X //,. (4.7) 

Using [«(x)|f = \u g (x') - i«W - " 3 (x')}| p < 2 p {|u 9 (x') - u(x)p + \u g (x')\P}, (4-7), 
and (4.6), we obtain the following: 



/ / x<s(x,x') \u g (x') - u(x)\ p dx'dx 

Bj Bj — 1 



>hj I Xs(x,x')\u(x)\ p dx'dx- f f X s(x,x')\ug(x')\P dx'dx 



2'' 



1 

2' J 



Xi(x,x') dx! \ \u(x) 

B j I J _ 1 j 



p dx- / J / X s(x,x')dx\ | U9 ( x ')| p dx' 

•/flj-i [jBj J 

= hB j _ 1 \[ |«(x)|* dx-l^-l / |^(x')| p dx' 
>^S m [ |«(x)|"dx- c out <5 m / 

z is,- -'B,_i 



a 

3 — ' 



Ug(x')\ P dx'. 



Hence, 



^ <5 m / |u(x)|f dx < a^(u g ,u) + c out S m / K(x')| p dx' 



2'' 

For the cases j = 1, 2, we respectively have: 



^ S m / |u(x)|" dx < o^K, u) + c ou t S m / | 5 (x')| p dx' (4.8) 
^ Jb x Jb q 

^S m f \u(x)\ p dx < a^(u g ,u) + c out S rn f \u{x')\ p dx' . (4.9) 

J Bo J- B\ 

In order to relate (4.9) to RHS of (4.8), multiply (4.8) by c ou t 

com 5 m I \u{x)\ p dx < C ^L^L a ^(u g ,u) + ^2lL S m j \g{x')\ p dx' . (4.10) 

J Bi ^in ^in J Bq 

Then, using (4.10), (4.9) becomes 



2p 



5 m / \u{x)\ p dx < a ip \u g ,u) + cout S m / \u{x')\ p dx' 



B2 J Bi 



<{^! + l} fl W( U9;!1 ) + {^ +Cra(}r / \g( X ')\ p dx'. 
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, I, using the fact that the covering of is 
-i-Bfc, Bj (~1 Bk = 0, j 7^ k, we arrive at the 



Repeating this procedure for j = 3, . . 
composed of disjoint strips, i.e., Q = LS l k 
result. □ 

Remark 4.2. For the condition number analysis, the underlying function space 
is chosen to be L2(£l). In addition, we enforce a zero Dirichlet boundary condition. 
Hence, the nonlocal Poincare inequality (4.2) reduces to the following: 



Ap„ CJ% 2 <5 m IM!l 2 (o) < a(u,u) 



(4.11) 



where a(u,u) := a^{u, u), and \p nC r.2 = <Vp racr , 2(m, fi) > with to > 1. 

Remark 4.3. When a general kernel function C is used, the Poincare constant 
depends on that function as well, i.e., \p ncr = Ap ncr (m,p, f2, C). This is the result 
given in [4, Prop. 2.5]. With the canonical kernel function in (4.1), the Poincare 
constant reduces to Xp n cr = Ap ncr (m,p, Q). Our main contribution is the incorpora- 
tion of the explicit 6 -characterization. This is the key step used in the analysis of the 
condition number and its dependence on the spatial dimension; see $4-2 and §^.5. 

4.2. The choice of Poincare parameter m and convergence to a local 
problem. After establishing a nonlocal Poincare inequality with 5-dependence explic- 
itly quantified, the big question becomes, What is to in (4.11)? The ^-quantification 
helps to answer that question, because we numerically observe that the value of m 
varies with the spatial dimension d; see Section 4.5. Next, we discuss a strategy to 
determine the value of to. This strategy is based on the fact that as 5 — > 0, a nonlocal 
bilinear form should converge to its corresponding local (classical) bilinear form. For 
discussions of convergence of other nonlocal operators to their classical local counter- 
parts, see [16, 36]. We resort to a Taylor series expansion after enforcing a sufficient 
regularity assumption. In fact, with u G C 4 (f2), we recover the numerically observed 
value of to = 3 for d = 1. We split the bilinear form in (3.5) into three pieces and 
introduce the change of variable in (3.6): 



a{u, u) 




u (x + e) de > dx + 5 u (x) dx 
S Jn 




u(x + e) de > u(x) dx 



=: T (1) +T (2) +T (3) . 
For the integrands, we use the following Taylor expansions: 



e du 



,2 j2 



d 2 u 



d 3 i 



u{x + £) = u{x)+ v.Tx + ».d^ + v.dx^ + oi£4) 



(4.12) 



2 e du 2e 2 



u (x + e) = u (x) + —u— + — <;( — ) +u — 



2! 



du 



dx 



dx 2 



2 e 3 f du d 2 u d 3 u 1 
^ 3 1- u. \ 



3! \ dx dx 2 dx 3 \ ^ ^ 



(4.13) 



Using (4.12) and (4.13), we obtain the following expressions for the integrals: 



T^ =6 I u 2 dx 




T (3) = _ 



28u 2 



fu 
dx 2 



2d 3 d 2 u 

^r u dx^ + oiS ' 



0{6 5 ) 



dx 



dx 



(4.14) 
(4.15) 
(4.16) 
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Combining (4.14), (4.15), and (4.16), denoting the local bilinear form by 



(u,u) := \u\% 1(n) , 



we obtain the following: 



a(u, u) 



8 3 ( du 



3 \dx 
£(u,u) + 0(S 5 ) 



0{6 5 ) 



dx 



Therefore, the scaled nonlocal bilinear form asymptotically converges to the local 
bilinear form: 



3 <T 3 a(u,u) = e(u,u) + 0(S 2 ). 
Using the nonlocal Poincarc inequality (4.11) and (4.17), we have 



(4.17) 



lim 3 A Pracr , 2 S m - 3 ||u||* <£(u,u). 

— >u 



Therefore, for the left hand side to remain finite, we have to enforce m > 3. We desire 
for the largest possible lower bound in the nonlocal Poincare inequality, and m = 3 
is the numerically observed value of m in ID; see the experiments in §4.5.1. 

4.3. An upper bound for a(u,u). Without the need of a discretization of 
a{u, u) as in the local case upper bound (i.e., inverse inequality), we prove the following 
dimension dependent estimate: 

Lemma 4.4. Let fl C M. d be bounded. Then, there exists A = A(O) > such that 



a(u,u)<XS d \\u\\l 2{ny 



(4.18) 



Proof. We provide the proof in 3D. The lower dimensional cases follow easily 
First, we utilize a uniform partition of Q composed of disjoint cubes of edge size 6. 
Namely, 



(li,liM 



n c 



u 



(j,j,fe)=(l,l,l) 



Hence, \b^j^)\ = <5 3 - For a fixed S, h,lj, and Ik are functions of £1 only. The proof 
relies on the important observation that 



{x' : llx-x'H <S, x G b (id>k) } C 



(i+l.j'+l.k+l) 

E 

(u,jj,kk)=(i-i,j-i,k-i) 



J (H,jj,kk) ■ 



(4.19) 



For ease of notation, we replace the canonical kernel x<5( x j x ') by 1 in its sup- 
port due to (4.1). Furthermore, we introduce the characteristic function \a to avoid 
integration outside of fl for the pieces of the cubes that may stay outside of f2. 
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We split the bilinear form in (3.5) into three pieces: 



a(u, u) = — f I ( l u ( x ') ~ u ( x )| 2 dx' > dx 

2 Jq I J\\x-x'\\<8 



< - I 



u(x')\ 2 dx'\ dx+- [ \u(x) 

I \ I H*!)\ dx'\\u{x)\ dx 

Jn J\\x-x'\\<8 



in I J||x-x'||<(5 
= : +I {2) +I (3) 



First, let 



r« 1 



^ m --=- 2 I / HxTdx' dx 

J ,b (i,j,k) I J l|x-x'||<<5 I 



• /b «,i,fe) l-'llx-x'||<a I 



lg) k) := / Kx')|dx' Kx)|dx. 



Then, 



/(l) + 7 (2) + j(3) = ( ' y ^ r(D + r(2) + r(3) 

(ij.fc) (i,3,k) (i,J,fc) 

(i,j,fc) = (l,l,l) 



We will give an upper bound for each piece. Using (4.19), we immediately see: 



(»+i.j+i,fc+i) 

r Sm^U E / Mx>(x')| 2 dx' dx 



< 



^ (ii,jj,kk) = (i-l,j-l,k-l) 
(i+l.j+l,fc+l) 

E llxn«lli,( 

(ii,jj,kk)=(i—l,j—l,k-l) 



Hence, we obtain an upper bound for JWby the following sum. 

i {1) = E < ^ imiW 
(t,j,fc)=(i,i,i) 
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Using (4.19), we immediately have: 

(i+l.j+i.k+i) 



I$ hk) <\i |xn(x).(x)| 2 dx £ f 

= o / lxn(x)w(x)| 2 dxi 
1 J hi,s*) { 



dx' 

W.i.fc) ^ (ii,jj,kk)=(i-l,j-l,k-l) " o (».«,*>*0 

(i+l.j+l,fc+l) 

E IW*>I 

(iijj,kk) = (i-l,j-l,k-l) 
(t+l.j+l,fc+l) 

= \\\xou\\^ k) \ E * 3 

^(M,ji,fcfc) = (j-l,j-l,fe-l) 

Hence, we obtain the upper bound for 1^ by the following sum. 

Ii2) - "if 43,^ c(2) * 3 HW 
(t,j,fc)=(i,i,i) 



Finally, using (4.19), Holder's inequality for each term and the identity 
< \{X 2 + Y 2 }, we obtain the following: 

(t+i.j+i,fc+i) 



< I IXn(x) u(x)| dx J I l^( x ') "( x ')l dx ' 

Jb li.3,k) [ (ii,jj,fefc) = (i-l,i-l,fc-l) Jb di,jj.kk) 

!(l+l.j + l,fc+l) 
E l 6 («^,fcfc)r /2 ||xn"||L 2 (& (ll ,„, fcfi) ) 
(ii,jj,kk)=(i— l,fc— 1) 

!(i+l.j + l,fc+l) 
E llxnulU,^. .,,) ||xnu|U 2 (6 (ji ^, fefe) ) 

(«,ji,fefc)=(t-l,3-l,fc-l) 
(t+l.j + l.fc+l) 



{ E Hxn u lli 1 i(W)) + llxn u lli 3 C6(«^,»*)) 

(ii,jj,kk)=(i— — l,fc— 1) 



Hence, we obtain an upper bound for i^ 3 - 1 by the following sum. 

i {3) = E < ° w ^ imiW 
(t,i,fc)=(i,i,i) 

The result follows from adding the upper bounds for the three pieces. □ 
Remark 4.5. The upper bound is numerically sharp when a piecewi.se constant 
discretization ofa(-,-) is used; see §^.<5. 

4.4. The Conditioning of the Stiffness Matrix K. Combining the nonlocal 
Poincare inequality (4.11) and the upper bound (4.18), we arrive at the condition 
number estimate: 

Theorem 4.1. The following spectral equivalence holds: 

\pncr.2 8 m <P^<\8 d . (4.20) 
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Let K be the stiffness matrix produced by the discretized a(u,u). Then, the condition 
number of K has the following bound: 

k(K) < 5 d - rn . (4.21) 



(a) Constant S, vary h. 



1/h 


S 


A mi „ # 


A ma x 


Condition # 


20 


0.3 


3.53E-03 


3.32E-02 


9.41E+00 


40 


0.3 


1.87E-03 


1.67E-02 


8.92E+00 


80 


0.3 


9.74E-04 


8.36E-03 


8.58E+00 


160 


0.3 


4.99E-04 


4.18E-03 


8.39E+00 


320 


0.3 


2.52E-04 


2.09E-03 


8.29E+00 


(b) Constant h, vary <5. 


1/h 


<5 


A mi „ 


A m ax 


Condition # 


100 


0.02 


2.90E-07 


4.50E-04 


1.55E+03 


100 


0.04 


2.13E-06 


9.54E-04 


4.49E+02 


100 


0.08 


1.66E-05 


1.93E-03 


1.16E+02 


100 


0.16 


1.29E-04 


3.80E-03 


2.93E+01 


100 


0.32 


9.36E-04 


7.03E-03 


7.51E+00 



Table 4.1 

Condition number for K in ID for (a) constant 8, allowing h to vary, and (b) constant h, 
allowing 8 to vary. This data is plotted in Figure J^.l. 




Fig. 4.1. Condition number for K in ID for (a) constants, allowing h to vary, and (b) constant 
h, allowing S to vary. The condition number is h-independent, but varies with S~ 2 . These figures 
are plotted from data in Table 4.1. 



4.5. Numerical Verification of Condition Number by a Finite Element 
Formulation. For all computational results in this paper, we assume that Q, = [0, l] d 
is the unit d-cube with a conforming triangulation 7^ where each element K of T~h is 
a d-cubc with a side length h > 0, where d is the spatial dimension. Consequently, 
each element in ID, 2D, and 3D is a line segment of length h, a square of area h 2 , 
and a cube of volume /i 3 , respectively. Let Vh C V be a finite dimensional subspacc 
of V from (3.1). We use a Galcrkin finite element formulation of (3.2): 



a(u h ,v h ) = (b h ,v h ) \/v h eV h . 
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(a) Constant 5, vary h. 



1/h 


s 


A mi „ # 


A n „ 


Condition # 


6 


0.333 


5.86E-04 


1.00E-02 


1.71E+01 


12 


0.333 


1.57E-04 


2.54E-03 


1.62E+01 


24 


0.333 


3.98E-05 


6.40E-04 


1.61E+01 


(b) Constant h, vary <5. 


1/h 


<5 


A mi „ # 


Amax 


Condition # 


24 


0.083 


3.22E-07 


4.05E-05 


1.26E+02 


24 


0.166 


4.32E-06 


1.67E-04 


3.86E+01 


24 


0.333 


3.98E-05 


6.40E-04 


1.61E+01 



Table 4.2 

Condition number for K in 2D for (a) constant 8, allowing h to vary, and (b) constant h, 
allowing 8 to vary. This data is plotted in Figure J i .l. 




Fig. 4.2. Condition number for K in 2D for (a) constants, allowing h to vary, and (b) constant 
h, allowing 5 to vary. The condition number is h-independent, but varies with 5 -3 / 2 . These figures 
are plotted from data in Table J h .l. 



For convenience, we have used Dirichlct boundary value for the theoretical construc- 
tion. To verify our theoretical results numerically, for implementation convenience, 
we utilize the discretization of a pure Neumann boundary value problem. 

a(u h ,v h ) = Vv h eV h , (4.22) 

where we take Vh to be the space of piccewise constant shape functions on the mesh 
T~h- We denote the resulting stiffness matrix by K. We numerically determine the 
largest and smallest nonzero eigenvalues, defining the effective condition number of 
the problem. 

4.5.1. Results in One Dimension. Results in this section appear in Tables 
4.1 and Figures 4.1. We first compute the condition number of K for different h 
while holding 5 constant to verify our conjecture that the condition number of K is 
/i-independcnt. The minimum nonzero and the maximum eigenvalues depend linearly 
on h, with a slope of unity. Consequently, the condition number of K is mesh size 
independent. We then compute the condition number of K for different <5 while 
holding h constant, and observe that the condition number varies with S~ 2 . Further, 
the maximum eigenvalue is proportional to S d (d = 1 in this case), in agreement 
with Lemma 4.4. Lastly, the minimum eigenvalue varies as 5 3 , in agreement with our 
finding of to = 3 in §4.2. This suggests that, in one dimension, we should redefine 
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r 



X 




BQ.2 



Fig. 5.1. Nonlocal two-subdomain problem. Note that the interface T is d = 2-dimensional. 



for consistency with the weak form of the classical (local) Laplace operator in the 
limit 8->0. 

It is instructive to compare convergence results for the nonlocal problem with the 
classical convergence results for the (local) discrete Laplace equation. For the discrete 
local Laplace equation, the condition number varies with h~ 2 [40, Theorem B.32]. In 
light of these numerical experiments, it appears that the characteristic length scale 8 
of the nonlocal problem (1.1) plays the role that the mesh size h plays in the discrete 
(local) Laplace equation. This is interesting, as S is a parameter of the continuum 
model, whereas ft. is a parameter of the discrete model. 

4.5.2. Results in Two Dimensions. Results in this section appear in Tables 
4.2 and Figures 4.2. We first compute the condition number of K for different h 
while holding 8 constant to verify our conjecture that the condition number of K is 
h- independent. The minimum and maximum eigenvalues depend linearly on h, with 
a slope of two, and again the condition number of K is mesh size independent. We 
then compute the condition number of K for different 8 while holding h constant, 
and observe that the condition number varies with S~ 3 / 2 . Further, the maximum 
eigenvalue is proportional to 8 d (d = 2 in this case), in agreement with Lemma 4.4. 

5. A Nonlocal Two-Domain Problem. We will construct both strong and 
weak (variational) formulations for domain decomposition. We first identify the pieces 
of the domain for this decomposition. Consider the domain in Figure 5.1. The nonlo- 
cal boundary of fi, BVL, is defined to be the closed region of thickness 8 surrounding Q. 
Let T be the open region corresponding to the interface between the two overlapping 
open subdomains f^ 1 ) and f^ 2 ). We define the overlapping subdomains i = 1,2, 

as the following: 



where 1^ is the open line segment adjacent to 0^ and T. Let Bfli be the nonlocal 
closed boundary of Qi that intersects Bfl. In addition, the remaining pieces Tn and 
Ts are the open regions to the north and the south of T, respectively. The nonlocal 
boundary of fiW is defined as follows: 



C(x, x') in (4.1) as 




otherwise. 




-.= fi,urur„ 
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The three equivalence results we will prove next are the main domain decompo- 
sition contributions of this article: 

1. Equivalence of the one-domain strong and two-domain strong forms. 

2. Equivalence of the one-domain weak and two-domain weak forms. 

3. Equivalence of the two-domain strong and two-domain weak forms. 

The fact that the main problem (1.1) is an integral equation provides a convenient 
framework for establishing these equivalences, especially for the last item. 

5.1. Two-Domain Strong Form. We seek a solution to (1.1) in L2(f2), i.e., 
u G Z,2(f2). We present a two-domain strong formulation of (1.1) and its equivalence 
to the original single-domain formulation (1.1). Let u^'(x) := Xqw (x) u(x). For the 
sake of brevity, we utilize the following notation: 

£n 4 (u)(x) := -xn 4 (x) / X*(x,x') • [u(x') - u(x)] dx' (5.1) 

C T (u)(x) := - Xr (x) f x*(x,x') • [u(x') - u(x)} dx'. (5.2) 

The two-domain version of (1.1) may be stated as follows: For % = 1,2, 

£n i (^ (i) )(x)=xn i (x) b(x), (5.3a) 
Xr (x) U W =Xr ( x ) U W, (5.3b) 

J] ^(x)/V(«W)(x)=^(x)6(x). (5.3c) 

»=1,2 

Remark 5.1. In /ig/if o/ (5.3), we make the following crucial observation about 
nonlocal domain decomposition. The strong version of the two-domain decomposition 
turns into a three-domain decomposition in which we produce three domain conditions, 
i.e., (5.3a) fori = 1,2, and (5.3c), together with the transmission condition (5.3b). 

We now show that the two-domain strong form is equivalent to the one-domain 
strong form. 

Lemma 5.2. The problems (1.1) and (5.3) are equivalent. 
Proof. (1.1) =*> (5.3) : 

Let u« = u| nw E V( l > , where is defined in (5.7). Using the nonoverlapping 
decomposition of the domain, we have: 

Xn(x) = Xfii(x) +xn 2 (x) + Xr( x )- ( 5 - 4 ) 
Using (5.4), suppressing the x-dependence, 

C(u) = xni £(u) + Xn 2 C{u) + - Xr ^i u ) + \ Xr 

- Xfl! + xn 2 £n 2 (« (2) ) + ^ /: r (" (1) ) + ^ ^r(« (2) ) = (5-5) 

Multiplying (5.5) by xn^ we g et (5.3a): 

xn, £n 4 (u w ) = xo ; 6. 
Likewise, multiplying (5.5) by Xf~; we S e t (5.3c): 

l -C T {u^)+ l -C T {u^) = Xrb. 
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Finally, we trivially obtain (5.3b) by the following: 



t (D|_ =u |_ = u (2)|_ 



(5.3) => (1.1) : 
Letting 



u := { U W, in fi 2 (5.6) 

{ uW -- 

and adding (5.3a) and (5.3c), we obtain: 




b = E fcni(« (i) ) + ^(« (i) ) 

i=l,2 ^ 

i=l,2 ^ 

= E I + 2 x r £ r ( u ) 

i=l,2 ^ 



□ 

5.2. Two-Domain Variational Form. We present a two-domain weak formu- 
lation of (3.2) and its equivalence to the original single-domain formulation (3.2). We 
define the spaces, i = 1,2, 

6L 3 (fiW) : «|boW =0}, (5.7) 
^°:={ B eI 2 (fi (i) ) : vbnwurur, = o} , 

A := {/i G ^2(r) : M = v \r for some suitable v S Z,2,o(fi)} , 
and the bilinear form <Xq(;)(u,v) : Z^^^) x Z,2(^ ) — > R as follows: 

o nW (u,u) := - / < / Xa~(x,x') • [u(x') - u(x)] dx' I • u(x) dx 

-\J T {j n ^ x ') • M x ') - W W1 dx '} • U W dx - ( 5 - 8 ) 

Remark 5.3. We have chosen the weights of ^ in the definition of the bilinear 
form in (5.8). Such a choice typically indicates that the underlying body is homoge- 
neous. For an inhomogeneous body, the weights can be chosen as a convex combination 
reflecting the inhomogeneity. 

We utilize the following notation to suppress the integrals in (5.8): 

a ni (u,v) := - |^x«(x,x')-Kx')-«(x)] dx'j • v(x) dx (5.9) 



10) 
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Wc can now represent the bilinear form (5.8) as: 

Cow (u, v) = i a r (u, v) + an z (u,v). 

Remark 5.4. The test function u, = v\n t S VW>°, i = 1,2 has its support only 
in Hi not fiW . Hence, we may reduce the bilinear form (5.8) to 

anw(u (i) ,Vt) =an i (w (i) ,Vi)- (5.11) 

Although, ar(u^ l \ Vi) may appear to create a coupling between the subdomains, no such 
coupling exists because Vi vanishes on T. Therefore, subdomain condition (5.12a) is 
an expression only for subdomain O™. 

Now, we state the two-domain weak form following the notation of [30]: Find 
u (i) e v®, i = 1,2: 

a Qli) (u^,v i ) = (b,v i ) ni V^eFW' , (5.12a) 

u {1) =u {2) onT, (5.12b) 

J2 a n(i) (u^,H^n) = (b,fi) r + Yl (&>ft (< V)n, e A. (5.12c) 

i=l,2 i=l,2 

where denotes any possible extension operator from r U T, to V^ 1 '. Next, wc 
will show that the one- and two-domain weak forms are equivalent. The proof for the 
local case can be found in [30, Lemma 1.2.1]. 

Lemma 5.5. The problems (3.2) and (5.12) are equivalent. 

Proof (3.2) =S> (5.12) : 
Let uW = tt|n(o € and «j = e yW>°, i = 1,2. Extend these functions by 
zero extension; 



««« := 



0, otherwise 

Vi, in 0, 

0, otherwise. 



By LHS of (3.2) and using Vi\r = 0: 

a(0«w« ,0iVi) = - f If X5(x,x') ■ [r^V^x') - 0«tt«(x)] dx' }- • f^(x) dx 



= aj2 I (u (l) ,w ! :) 



By RHS of (3.2), 

(b,8iVi) = (MtW 

Hence, (5.12a) is satisfied. (5.12b) is trivially satisfied. 
Further, for fi 6 A define the function IZfj. as: 



Tlji := 



18 



Burak Aksoylu and Michael L. Parks 



Since IZfi lives only in U Fx U L U T 2 U 17 2 , it vanishes on SO. Therefore, lZ[i G V. 

From (3.2), partitioning the outer integral and using = n on r, 

we obtain the LHS of (5.12c): 

a(uM = \ a T {u^\ix) + \ a v (vP-\p) + ^ an, (u {i \ 7^V) 

i=l,2 

= i ar ( u a>,-K<iV) + I ar (u< a >,rc< 2 V) + X] «^(« W ,^ ; V) 

i=l,2 

Likewise, from (3.2) and partitioning the integral, we obtain the RHS of (5.12c): 

(b,Hii) a = (ft.^Vjti! + (&,ft (2 V)n 2 + (6,A*)r- 
Hence, we obtain the transmission condition (5.12c). 

(5.12) => (3.2) : 

Let Mr := (due to (5.12b), we also have ur = u^) and 

( infix 
tt := < u< 2 ), in fi 2 (5.13) 
[ Mr, in r. 

We partition the outer integral, use (5.13) and the transmission condition (5.12b). 
Then, for v G V, LHS in (3.2) becomes the following: 

i=l,2 

= ^ ar("r,v) + ^ a r (ur, «) + 5Z an<(u^,u) 

i=l,2 

= - a r (u (1) ,iO + - a r (u (2) ,w) + 53 a!),(M (t) ,D) 

i=l,2 

= 53 a nW («««). (5.14) 

i=l,2 

Let n := v\ T . Then, v-TZ^fx G yM<°. First, we add and subtract fj, to (5.14) and 
apply the domain conditions (5.12a) for v — TZ^fi. Then, we apply the transmission 
condition (5.12c) and use v\r = £t. Hence, we arrive at the RHS in (3.2): 

J2 a nw (u^,v) = a n w(« (i) ,i>-ft (< V) + 53 «qw(u W ,^V) 

i=l,2 i=l,2 t=l,2 

= 53(6, w -^) M ) fii + 5] a nW (nW,^/x) 

i=l,2 i=l,2 

= 53 (6, v - ■RPtin, + ( b > K®l*)ni + ( b > M)r 

i=l,2 i=l,2 

= 5Z + (M)r 

i=l,2 

= (M). 
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5.3. Equivalence of Two-Domain Strong and Two-Domain Weak Forms. 

Wc prove the equivalence of the strong and weak forms of the two-domain formulation. 

Lemma 5.6. The problems (5.12) and (5.3) are equivalent. 

Proof. (5.3) =*> (5.12) : 
Testing (5.3a) against u, £ V^'° and using = 0, wc obtain (5.12a): 

(Mi)fii = (Xfii £n<(« W ),«i) 
= «aw (w W , «,:)• 

(5.3b) and (5.12b) are identical. 

Next we show how to obtain (5.12c). Given fx £ A. Testing (5.3a) against TZ^'fi, 
wc immediately obtain: 

an 4 (u W ,^ W A*) = (b,TZ^n) ni . (5.15) 
Likewise, testing (5.3c) against we get: 

l - (x r ^ (1) ),M) + ^ (xr^ (2) ),f) = (Xr^)=^M)r (5-16) 
bmce is any extension operator, fi\r = H, we can rewrite (5.16) as follows: 

\ (x r ^ (1) ),^ (1) M) + 5 (Xr^ (2) )^ (2 V) = (&,M)r- (5-17) 
Adding (5.15) and (5.17), we obtain (5.12c). 
(5.12) =>■ (5.3) : 

Since for Vi £ V^'°, we have Uj|p = 0, (5.12a) reduces to the following: 

(Mi)n< = a S!,( uW .fi) 

= (£o I (^ (4) ),^). 

Equivalently, we have 

We obtain (5.3a). 

Since we have (5.3a), test it against TZ^'fx: 

(£n i (« (i) ),^ (i V) = (xn i &,^ (?; V) 
an i (u (i) ,^ i V) = (&,^ (i V)n i - (5-18) 

Subtracting (5.18) from (5.12c), we obtain 

\ a r (u^\n^fi) + l - aM 2 \n^ii) = (& )M )r. 
Using TZ^ /i\r = \i for the above equation, we arrive at: 

(i xr ^(« (1) ) + 5 XT ^ {2) ) - Xrb, m) = V,i£ A. 
We obtain (5.3c). □ 
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6. Towards Nonlocal Substructuring. Here we write out the linear algebraic 
representations arising from the two-subdomain weak form (5.12), identifying the dis- 
crete subdomain equations and transmission conditions. We then construct a nonlocal 
Schur complement, discuss its condition number as a function of h, 5, and provide 
supporting numerical experiments. 

6.1. Linear Algebraic Representations. We consider a finite element dis- 
cretization of (5.12). Letting vS % ' denote the finite element space corresponding to 
, we define: 



V, 



ZVf : ^Isn.urur, =0} 



{v h - *' W 

Ah := {fin ■ fJ-h = v h\r for some suitable v h G V h ] ■ 
We see that the finite element formulation of (5.12) can be written as: 



(1) (2) 
U h = U h 



Vui.h G V, 



(0,o 



i=l,2 



i=l,2 



„ h denotes any possible extension operator from U T 



on r, 
V/ih G A/j. 

K0 



i t h to V,) l; . Following 



where 

standard practice, we take these extension operators to be the finite element inter- 
polant, which is defined to be equal to [ih at the nodes in the thick interface T and zero 
on the internal nodes of fl;. If we number nodes in Oi first, nodes in Q 2 second, and 
nodes in T last, we will arrive at a global stiffness matrix that takes the traditional 
block arrowhead form: 



(6.2) 



6.2. Energy Minimizing Extension and the Schur Complement Condi- 
tioning. In order to study the conditioning of the Schur complement in the nonlocal 
setting, we first define an analog of the harmonic extension in the local case. 

Definition 6.1. For a given q E L 2 (T), E l (q) : L 2 (T) — > L 2>0 (Cl^) is called its 
energy minimizing extension into 0,, if the following holds: 





' K u 





K ir ' 




U\ 




' fi ' 


K = 







K 2T 




u 2 




h 




K rx 


K T2 


K rr 




u T 




fr 



Ei(q)\i 



q, 



(6.3) 



and a(Ei(q), v) = for all v G L 2i o(Q^') such that v\r = 0. 

The energy minimizing extension Ei(q) of q defines a canonical bilinear form 
Si(q,q) : L 2 (T) x L2(r) — > R that is associated to the interface T whose discretization 
corresponds to the subdomain Schur complement matrix below. 



Si(q,q) := a(E l (q),E. l (q)) 
(S®q h ,q h ) :=a{E % {q h ),E. l {qh))- 



(6.4) 
(6.5) 



Let us denote the restriction of u G L 2 ${Vl^) to Y by ur '■= u\r- The following 
discussion will reveal the reason why Ei{ur) is called an energy minimizing extension. 
Let us consider the following decomposition of u: 



u=[u- Ei(u r )] + E t (ur)- 



(6.6) 
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Since [u — E^uy)) |r = 0, by Definition 6.1 we have: 

a{u-E i (u r ),E i (u T )) = 0. (6.7) 

Using (6.6) and (6.7), we have the energy minimizing property of Ei(ur) among 
u £ -^2,0 (^ ) with u\r = ur'. 

a(u,u) = a(u - Ei(u T ),u- E t (u r )) + 2a(u - E t (u r ), E^up)) + a(E t (u r ), Ei(u T )) 
>a(Ei(u r ),Ei(u r )). (6.8) 

Therefore, using (6.8), (6.4) and (4.18), we have: 

s t (u r , u r ) < a(u, u) <XS d ||ti||i 3(n ci))) 

for all u G Li2,o(& l ')i m particular, for u xr- Hence, 

Sl (ur,u r ) < A 6 d ||u||£ a(r) . (6.9) 

For the lower bound, we simply use (6.3) and (4.11): 

Ap„ cr ,2<HMIi 2( r) < Ap„ CJ ..2(5 m ||^(ur)|li 2( c2) < a(E(u r ), E(u r )) = s{u r ,u r ). 

(6.10) 

We have proved the following spectral equivalence result: 
Theorem 6.1. For any q e L 2 (T), we have: 

Apncr.2 S m < * f 2 ' q) <X5 d . (6.11) 

Thus, the condition number of the Schur complement matrix St '■= S^ + has 
the following bound: 

<S r ) < S d - rn . 

Remark 6.1. The preceeding condition number estimate indicates that the con- 
dition number of the Schur complement is no greater than that of the corresponding 
stiffness matrix; see (4.21). This estimate is not tight. In fact, we numerically observe 
smaller condition numbers for the Schur complement; see Table 6.3. 

6.2.1. The Nonlocal Schur Complement Matrix. When the contributions 
from each subdomain are accounted separately, we can write Krr in (6.2) as Krr = 
JC« +K$. Then, £T« in (6.5) can be written as follows: 

S« := Afr - KnK^Kv. 

The solution across the whole of T is determined by solving SrUr = / for Ur, where 

/ := f r - - K T2 K^} 2 . 

We proved in §4.4 and verified numerically in §4.5 that the condition number of 
the stiffness matrix K is /i-independent. Therefore the condition number of the Schur 
complement matrix Sr should also be /i-indepcndcnt. We will examine this conjecture 
numerically in §6.3. 
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(a) Constant 8, vary h. 



1/h 


s 


Vin # 


A m „ 


Condition # 


20 


0.3 


2.18E-02 


3.05E-02 


1.40E+00 


40 


0.3 


1.15E-02 


1.52E-02 


1.32E+00 


80 


0.3 


5.91E-03 


7.55E-03 


1.28E+00 


160 


0.3 


2.99E-03 


3.77E-03 


1.26E+00 


320 


0.3 


1.51E-03 


1.88E-03 


1.25E+00 


(b) Constant h, vary S. 


1/h 


6 


Amin ^ 


A m ax 


Condition # 


100 


0.02 


2.37E-04 


4.10E-04 


1.73E+00 


100 


0.04 


5.49E-04 


8.17E-04 


1.49E+00 


100 


0.08 


1.20E-03 


1.63E-03 


1.36E+00 


100 


0.16 


2.49E-03 


3.23E-03 


1.30E+00 


100 


0.32 


5.07E-03 


6.44E-03 


1.27E+00 



Table 6.1 

Condition number for Sr in ID for (a) constant 8, allowing h to vary, and (b) constant h, 
allowing 5 to vary. This data is plotted in Figure 6.1. 




Fig. 6.1. Condition number for Sr in ID for (a) constant 8, allowing h to vary, and (b) 
constant h, allowing 8 to vary. The condition number of Sr is apparently nearly independent of 
both h and 8 in ID. These figures are plotted from data in Table 6.1. 



6.3. Numerical Verification of the Schur Complement Conditioning. To 

test the conjecture of the previous section, we discretize the pure Neumann boundary 
value problem 

s(u,u) = o xer, (6.12) 

where d is the spatial dimension, using piecewisc constant shape functions on uniform 
cartesian mesh, and numerically determine the largest and smallest nonzero eigenval- 
ues, defining the effective condition number of the problem. 

6.3.1. Results in One Dimension. We define the regions fJi = (0, 0.5 — 8/2), 
n 2 = (0.5 + 5/2, 1), and T = (0.5 - (5/2,0.5 + 8/2), such that T is always a region 
of width (5 centered at x = 0.5. We then compute the largest and smallest nonzero 
eigenvalues of Sr ■ 

We first compute the condition number of Sr for different h while holding 8 
constant to verify our conjecture that the condition number of Sr is /i-independcnt. 
Our results appear in Tables 6.1 and Figures 6.1. The minimum and maximum 
eigenvalues depend linearly on h, with a slope of unity. Consequently, the condition 
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(a) Constant 8, vary h. 



l/h 


s 


A mi „ # 


Vox 


Condition # 


6 


0.333 


1.34609E-03 


9.96E-03 


7.40E+00 


12 


0.333 


3.88334E-04 


2.52E-03 


6.48E+00 


24 


0.333 


1.09824E-04 


6.32E-04 


5.75E+00 



(b) Constant h, vary <5. 



l/h 


<5 


A mi „ # 


A n „ 


Condition # 


24 


0.083 


1.70E-06 


4.01E-05 


2.36E+01 


24 


0.167 


1.44E-05 


1.65E-04 


1.14E+01 


24 


0.333 


1.10E-04 


6.32E-04 


5.75E+00 



Table 6.2 

Condition number for Sr in 2D for (a) constant 8, allowing h to vary, and (b) constant h, 
allowing 8 to vary. This data is plotted in Figure 6.2. 




Fig. 6.2. Condition number for Sr in 2D for (a) constant 8, allowing h to vary, and (b) 
constant h, allowing 8 to vary. The condition number of Sr in 2D is nearly independent of h but 
varies with 5 _1 . These figures are plotted from data in Table 6.2. 



(a) ID. 



(b) 2D. 



l/h 


S 


Cond. # K 


Cond. # S r 




l/h 


5 


Cond. # K 


Cond.# S r 


20 


0.30 


9.41E+00 


1.40E+00 




6 


0.333 


1.71E+01 


7.40E+00 


40 


0.30 


8.92E+00 


1.32E+00 




12 


0.333 


1.62E+01 


6.48E+00 


80 


0.30 


8.58E+00 


1.28E+00 




24 


0.333 


1.61E+01 


5.75E+00 


160 


0.30 


8.39E+00 


1.26E+00 




24 


0.083 


1.26E+02 


2.36E+01 


320 


0.30 


8.29E+00 


1.25E+00 




24 


0.167 


3.86E+01 


1.14E+01 


100 


0.02 


1.55E+03 


1.73E+00 




24 


0.333 


1.61E+01 


5.75E+00 


100 


0.04 


4.49E+02 


1.49E+00 










100 


0.08 


1.16E+02 


1.36E+00 










100 


0.16 


2.93E+01 


1.30E+00 










100 


0.32 


7.51E+00 


1.27E+00 











Table 6.3 

Condition numbers of K and Sr in ID and 2D. The condition number of the Schur complement 
is less than the condition number of the original system. This data is a summary of data from 
Tables 4.1, 4.2, 6.1, and 6.2. 



number of Sr is /i-indcpcndcnt. We then compute the condition number of Sr for 
different 5 while holding h constant, and observe that the condition number is nearly 
<5 independent. 

The condition number of the Schur complement of the weak classical (local) 
Laplace operator varies as (Hh) -1 , where H is the subdomain size [40, Lemma 4.1.1]. 
Assuming S "plays the role" of h as in §4.5.1, one might expect that the condition 
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number of the Schur complement for the weak nonlocal Laplace operator to vary as 
(_ff<5) _1 , but Figure 6.1(b) shows that this is clearly not the case in ID. 

6.3.2. Results in Two Dimensions. We define the regions = (0, 0.5 — 
5/2) x (0, 1), Q 2 = (0.5 + 5/2, 1) x (0, 1), and V = (0.5 - 5/2, 0.5 + 5/2) x (0, 1), such 
that r is always a region of width 5 centered at x = 0.5. We then compute the largest 
and smallest nonzero eigenvalues of Sr- 

We first compute the condition number of Sr for different h while holding 5 
constant to verify our conjecture that the condition number of Sr is /i-indcpcndent. 
Our results appear in Tables 6.2 and Figures 6.2. The minimum and maximum 
eigenvalues depend linearly on h, with a slope of approximately two. Consequently, 
the condition number of Sr is nearly h- independent. We then compute the condition 
number of Sr for different 5 while holding h constant, and observe that the condition 
number varies as <5 _1 . 



Dim 


A" 1 (K) 

mm v ' 




k{K) 




AmaxCSr) 


«(Sr) 


ID 


o{s- [i ) 


0(S) 




ois- 1 ) 


0{8) 


0(1) 


2D 


0(8- 3 - 5 ) 


0(S 2 ) 


ois- 1 - 5 ) 


0(5- 3 ) 


0(S 2 ) 


Op" 1 ) 



Table 7.1 

The 8 -quantification of the reported numerical results. 



7. Conclusions and Future Work. We collect the numerical results in Ta- 
ble 7.1. We observe that k(K) is independent of the mesh size, and that and n(Sr) 
is nearly so. The numerical estimates are computed using a piecewise constant dis- 
cretization. Therefore, in order to show the sharpness of the estimates analytically, 
one can use the same piecewise constant functions. This is a subject for future work. 

We conclude the sharpness of the upper bound 0(5 d ) in the spectral equivalence 
for a(u,u) and s(u,u), i.e., (4.20) and (6.11), respectively. When the condition num- 
bers of the nonlocal and local problems are compared, we see that k(K) = 0(5~ 2 ) 
is analogous to the local case, where the condition number varies as Oihr 2 ). How- 
ever, this does not hold for higher dimensions because k(K) = 0(5 d ~ m ), and hence, 
depends on the spatial dimension unlike in the local case. By demonstrating the non- 
local problem converges to the local one as 5 — > 0, we showed that m = 3 is sharp in 
ID. Numerically we observed that m ~ 3.5 in 2D. Also, the fractional power is some- 
what surprising. A systematic and constructive mechanism for obtaining the values 
of m = m{d) in higher dimensions is unclear. Therefore, more comprehensive study 
is needed to explore the dependence of m upon, for instance, the underlying basis 
functions and spatial dimension. Such study is also needed for the Schur complement 
to obtain sharp bounds. The effect of subdomain sizes on the conditioning of the 
Schur complement is unknown. 

Application of an appropriate preconditioner, involving the solution of a coarse 
problem, reduces the condition number of the Schur complement of the weak classical 
(local) Laplace operator from O ((Hh)- 1 ) to O ((1 + \og{H/h)) 2 ) [40, Lemma 4.11], 
[39, §4.3.6]. One unexplored area involves examining the role of a coarse problem 
in the nonlocal setting, which has not been considered here. A logical direction 
would be to expand other substructuring methods to a nonlocal setting, such as 
Ncumann-Dirichlct, Neumann-Neumann, FETI-DP (the dual-primal finite element 
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tearing and interconnecting method) [26], or BDDC (balancing domain decomposition 
by constraints) [15]. Additional opportunities for future research include addressing 
convergence analysis for alternative domain decomposition methods not based on 
substructuring in a nonlocal setting. More fundamental concepts in Schwarz theory 
such as stable decompositions and local solvers need to be reconstructed for nonlocal 
problems to support convergence analysis for additive, multiplicative, and hybrid 
algorithms. 

Appendix A. The Poincare constant is independent of S. 

In order to show that the Poincare constant Xp ncr is independent of 5, we give 
the construction in detail that leads to (4.6). Let the radius of inscribing and circum- 
scribing spheres of ft be denoted Ti n and r out , respectively. We compare the volume 
of Bj to the volume of the annuli of width of the inscribing and circumscribing 
spheres. 

In 2D, we give the explicit expression of the area of circumscribing annulus and, 
for sufficiently small S, we have the following: 
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In 3D, the volume of the inscribing annulus is given by the following: 
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Similarly, \B 
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Hence, 



Crnrf-H" 1 < |% m | < \Bj\ < \B jiOUt \ <c^r, 



d-l 



(A.l) 



Note that both c in rf„ 1 
dimension, we conclude (4.6), i.e., a n S r 
only on Q. 



and c out r d ou t depend only on J7. Consequently, in any spatial 
< \Bj \ < c ou t 8 m , where Ci n and c ou t depend 
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